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We demonstrate a tunneling and rectification behavior in bilayer graphene. A 
stepped dielectric top gate creates a spatially modulated electric held, which opens 
the band gap in the graphene and produces an insulating region at the p-n interface. 
A current-voltage relationship exhibiting differential resistance peak at forward bias 
stems from the tunneling current through the insulating region at the p-n interface. 
The tunneling current reflects singularities in the density of states modified by the 
electric field. This work suggests that the effect of carrier charge tuning by external 
electric field in 2D semiconductors is analogously to that by impurity doping in 3D 
semiconductors. 
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I. INTRODUCTION 

In an atomically-thin two-dimensional (2D) conductor like graphene, carrier charges can 
be induced with an external electric field PQ. Inducing carrier charges by an electric field 
is one significant advantage of 2D conductors, because it provides uniformly distributed 
carrier charges and overcomes the problem of random spatial fluctuation of dopants|2j. This 
random fluctuation is one of the important problems in nanoelectronics that use traditional 
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intrinsic semiconductors. A top gate that induces distinct localized electric fields in the 2D 
conductor can produce distinct regions with different carriers, thus forming a p-n junction 
in graphene that are required for devices[2l 0]. Klein tunneling^ E] arid quantum Hall 
edge modesjTJ |8] have been observed in these junctions. In monolayer graphene, such gate- 
controlled p-n junctions have been demonstrated increase the resistance through the active 
regions. However, an efficient barrier to current flow that will produce diode current-voltage 
characteristics has not yet been reported. An efficient barrier requires a potential barrier 
region to be induced in a 2D conductor with a bandgap. 

In this paper, we report the realization of such a barrier region at a p-n junction in bilayer 
graphene (BLG) , the first observation of tunneling through a p-n junction induced by electric 
fields in semiconducting BLG, and the first instance of rectification in a device based on a 
2D material. We observed a differential resistance peak at forward source-drain bias. We 
attribute this peak to tunneling between p and n regions, and credit it for diode behavior 
similar to what is observed in a 3D Esaki diode [SHU]. We employed a uniform bottom gate 
and a stepwise top gate to form a p-n interface. A thin top gate dielectric (less than 10 nm) 
[12] is critical to shrink the spatial transition between p and n regions. The electric field 
between the top gate and the bottom gate opens a band gap in the BLG [T3TH8] . making a 
tunnel barrier between p and n regions. As this is the first instance of rectification in 2D, we 
anticipate our approach will provide a starting point for creating gate-controlled diodes in 2D 
conductors. Furthermore, a gate-controlled p-n junction would be utilized in optoelectronic 
devices operating in the THz regime [121 [20J, which is covered by the gate-tunable range 
of the band gap in BLG. Utilizing the tunneling effect has an advantage for high speed 
electronics which would be one major application area for graphene electronics[21J. 



II. EXPERIMENTAL APPROACH 



An electric field produced by a stepped-gate can generate opposing charges under each 
half of the gate (Fig. 1(a)). A charge-neutral region exists under the step. If the 2D 
semiconductor has a band gap, charge carriers are depleted in the charge-neutral region. 
The depletion region has an in-plane electric field which originates from the gate electric 
field, canceled partially by the electric field from charge carriers in the 2D semiconductor. 
This mechanism sharply contrasts the impurity- doped p-n junction in a 3D semiconductor, 
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where carrier recombination produces a depletion region. The depletion of carriers leaves 
charged donor and acceptor impurities, resulting in a built-in electric fieldjD]. We note 
that the charge-neutral region and charged region in the gate-controlled p-n junction is 
opposite from the impurity-doped one: In the gate-controlled one, the depletion region is 
charge-neutral, while p and n regions are charged by carriers. In the impurity-doped one, 
the depletion region is charged by dopants (donors and acceptors), and uniform regions are 
charge-neutral. In spite of these differences, a model based on electrostatics suggests that 
the gate-controlled p-n junction in 2D mirrors the operation of the impurity-doped one in 
3D (Appendix A). 



III. SAMPLE FABRICATION 



BLG samples were prepared from kish graphite by the mechanical cleavage method[T] and 
adhered on to a highly-doped Si substrate with a 90-nm-thick SiC>2 surface layer. A graphene 
sample with multiple electrodes was patterned by oxygen plasma etching for four-terminal 
measurements (Fig. 1(b)). The bilayer channel was 0.4 fim in width and was sandwiched 
between a substrate bottom gate and a stepped top gate. The top gate was composed of two 
regions with different gate dielectric thicknesses. Half of the area of the graphene channel 
(surrounded by dashed lines in Fig. 1(b)) was covered by a 5-nm-thick layer of SiCV Then, 
the entire area of the graphene between the voltage terminals was covered by a 30-nm-thick 
Al film. The sample was exposed to air for several hours for partial oxidization of the Al 
film. An oxidized (A10 x ) layer formed not only on the surface but also at the interfaces of 
Al/ graphene [TBI 1221 [23] and Al/Si0 2 [2U [25]. The Si0 2 layer between part of the graphene 
and the A\O x layer, increases the dielectric thickness over that region [21] which results in 
the formation of a stepwise junction in the top gate (Fig. 1(c)). When a voltage is applied 
to the top gate, two different electric fields are simultaneously applied to the individual 
regions. The field effect mobilities extracted from the gate voltage dependence were 1300 
cm~ 2 V _1 sec _1 for electrons and 1800 cm~ 2 V _1 sec _1 for holes, independent of the dielectric 
thickness. 
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IV. RESULTS AND DISCUSSION 

The two gate regions with different dielectrics thicknesses create charge neutrality point 
(CNP) ridges with different slopes in the Hg-K g plane. Figure 1(d) shows a contour plot of 
the resistance Ro through the junction as a function of Vbg and V tg . The two ridges separate 
the Vbg-Vtg plane into p-p, p-n, n-p, or n-n combinations of carrier polarity in the graphene. 
On the ridge of CNPs, the height of the resistance peak increased with an increasingly 
negative electric field (Vb g —> —40 V) as an evidence of the band gap opening under the 
electric field [IB]. From the slope of the CNP lines, we can extract the top gate capacitance 
values [3J: C tg i = 9.3 x 10" 3 F/m 2 in the thick region and C tg2 = 5.5 x 1CT 3 F/m 2 in the thin 
region. If we assume the dielectric film is simply composed of SiC>2, equivalent thicknesses 
are approximately 3.7 nm (= di) and 6.3 nm (= c? 2 ), respectively. The stepwise top gate 
causes a transient region of ~ 5 nm which gives a width of the insulating region Wd (Fig. 
1(e) and Appendix A). 

In an actual bilayer of graphene with a band gap, a disordered potential forms a band 
tail [27]. The band structure forms two peaks in the density of states (DOS) at the edge of 
the valence band and the conduction band[2El [29]. The energy gap between the two peaks 
.Eg hardly depends on disorder, and is almost identical to the band gap in the unperturbed 
system [3D]. This gap can be extracted from the temperature dependence of the conductance 
at the CNP[18j. The CNP conductance comprises the intrinsic band conduction and the 
hopping conduction via localized states which make up the band tail. The former has a 
thermal activation energy E g /2. Extracted E g is ~ 0.2 eV at Vb g = —40 V both in the thin 
and thick dielectric region. The band gap makes an insulating region inserted between the p 
and n regions. Using this band gap, we estimate tunneling probability across the junction to 
be 0.07 and the tunnel resistance R t to be several kfl (Appendix B). The hopping conduction 
causes leakage current coexisting with the tunnel conduction. The leakage resistance is 
estimated to be ~ 1 kQ (Appendix B), which is comparable to the tunnel resistance. 
Then, the junction resistance, R s = 1/(1/R t + 1/Ri), is on the order of 1 kO. Using this 
value, it is estimated that the voltage drop at the junction Vj is a few percent of Vd s - 

The differential resistance dV^/dl was measured in a four-terminal configuration across 
the p-n junction as a function of the source-drain voltage Vd s , rather than the leakage 
current, to investigate the junction property because the differential resistance is sensitive 
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to a nonlinear tunneling current. We found a dVd s /dI peak for a large bottom gate voltage 
(Vbg = — 40 V) that opens the band gap (Fig. 2(a)). The peak was not observed for a bottom 
gate voltage (Vb g = —32 V) that was too small to open the band gap (bottom-right inset of 
Fig. 2(a)). The increase of the dVd s /dI at large Vd s is caused by charge redistribution by the 
Vds bias[3Tj. The peak appeared at a forward bias of Vd s ~ 50 mV, regardless of Vb g or V tg 
(Fig. 2(b)). The peak height depends on the bottom gate voltage and became pronounced 
when the gate electric field was increased by applying V^ g (Fig. 2(c)). 

We analyze the dV^s/dl peak observed in the experiment. A DOS diagram of the un- 
biased p-n junction is illustrated in Fig. 3(a)i). The voltage drop, Vj, depresses the 
energy in the p-type side (Fig. 3(a)ii)). The tunneling current, I t , from the p to n 
region is given by J t oc T t / [f(E - E F ) - f(E - E F + eVj)] D n (E)D p (E + eV^dE, where 
f(x) = 1/ [exp(x//cB^) + 1] is the Fermi distribution function, and D p (e) and D n (e) are the 
DOS at the energy level e (measured from the mid-gap) for the p-type and n-type sides [5]. 
The tunneling current reaches a maximum when the peaks in D p (E + eVj) and D n (E) align 
at the energy level (Fig. 3(a)ii)). For a larger Vj, the tunneling current becomes smaller be- 
cause the DOS peaks go out of alignment (Fig. 3(a)iii)). As a result, the tunneling current 
has a peak that is a function of Vd s - In the total current, the tunneling current becomes 
indistinguishable from the leakage current via localized states when the leakage current is 
comparable to or larger than the tunneling current. The peak structure displayed a dV& s /dI 
peak, as shown in Fig. 3(b), which was observed at Vd s ~ 50 mV. Because the voltage drop 
at the junction is a few percent of Vd s , the energy difference between the Fermi level and 
the band edge must be on the order of a few meV. 

The temperature dependence of the dV& s /dI peak (Fig. 4(a)) gives reasonable sup- 
port to the model of the tunnel junction discussed above. Using the Sommer- 
feld expansion in terms of temperature [32j, the tunneling current is proportional to 
T t JEp-eVj + a i(^B^) 2 -P"(-^)] dE within the second order of the temperature T, where 

P(E) = D n (E)D p (E + eVj), P"(E) = d 2 P(E)/dE 2 , and a x = \ JZo x 2 {-£^y^) dx ~ 
1.6. The function P(E) represents a density of states for elastic tunneling between the p 
and n region. Because a temperature coefficient for the second order is proportional to 
the integration of P"(E) around the Fermi level, a trend of the tunneling current in the 
temperature dependence is determined by convex upward or downward in the P(E). For a 
peak in tunneling current, P(E) has a peak around the Fermi level, i.e. P"(Ep) < (Fig. 
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4(b)). Thus, the tunneling current decreases with increasing T. This trend can be observed 
in the temperature dependence of dVd s /dI-Vd s , as shown in Fig. 4(a); the decrease in the 
tunneling current peak causes the dV^/dl peak to diminish as the temperature increases 
(Fig. 4(c)). The decrease in the tunneling current peak also leads to the Vd s dependence 
in the dV& s /dI-T curve (Fig. 4(d)). At a low Vd s bias (Vd s = 0), dV^s/dl becomes larger at 
higher temperatures, reflecting the decrease in the tunneling current. In contrast, at higher 
Vd s bias voltages, dV& s /dI becomes smaller at higher temperatures, reflecting the decrease 
in the dV& s /dI peak (Fig. 4(c)). Here, a second-order temperature coefficient (1.3 x 10 -5 
K~ 2 ) is extracted from the dV^/dl-T curve at the dV^/dl maximum (Vd s = 50 mV). A 
corresponding energy for the temperature coefficient is 20 meV. Because the temperature 
coefficient (oc P"(Ep)) determines the sharpness of the P{E) peak, the corresponding energy 
represents a broadness of the peak, caused by band tail of the DOS. This indicates that the 
band tail width is approximately 20 meV which on the order of 10% of the energy gap E g . 

The thermal energy at higher temperatures also generates a leakage current via localized 
states. The leakage current becomes larger with higher T, while the tunneling current 
becomes smaller. The tunneling current is more dominant than the leakage current in 
the temperature dependence, because the dVd s /dI at Vd s = becomes larger at higher 
temperature. The overall trend of the dV<\ s /dI peak observed in this experiment is similar 
to the negative differential resistance (NDR) of the Esaki diode [9HTT]. For a larger NDR 
in graphene, the leakage current must be suppressed. The impurity states in the band gap 
under high electric field are the most probable cause of the leakage current. Thus, reducing 
impurities and defects is important for fabricating a diode device governed only by the 
tunneling effect. 

V. CONCLUSIONS 

We observed clear tunneling signals in semiconducting BLG p-n junctions. This provides 
the first experimental evidence for the existence of conduction barrier in gap tunable atomic- 
layer conductors, which is essential to realize wavelength-tunable optoelectronic devices. 
We also identified localized states as the source of the diffusion current within the band 
gap, which highlights the importance to exclude the impurities and disorder in graphene to 
improve performance. With appropriate biasing conditions and transparent top gate stacks 
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for THz electromagnetic waves, novel wavelength-tunable optoelectronic devices would be 
viable in BLG. 
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Appendix A: Width of the Depletion Region 

We derive a Poisson equation to describe the electric potential profile in a double-gated 2D 
semiconductor. Here, we define a 2D conductor as a thin conductor which is much thinner 
than electric field screening length. Monolayer and bilayer graphene fit to the definition [22J. 
In the 2D semiconductor, the charge distribution along the thickness direction is negligible. 
The 2D semiconductor (thickness d s , dielectric constant e s ) is sandwiched between a bottom 
gate and a top gate (Fig. A- 1(a)). We consider a small region in the semiconductor from 
a position x to a position x + Ax. Gauss's law gives a relationship among the top gate 
electric field (E t ), the bottom gate electric field (E^), the electric field in the semiconductor 
plane (£■;), and charge density in the small region q(x): e s (Ei(x + Ax) — E{(x)) d s + (s t E t — 
e^E h )Ax = q(x)Ax, with e t and as the dielectric constant of the top gate insulator and 
the bottom gate insulator, respectively. Taking the limit as Ax approaches zero, we have a 
one-dimensional Poisson equation for the electric potential ( &) in the 2D semiconductor, 

d 2 ^ _ _dEi _ e t E t - e h E h - q(x) 
dx 2 dx e s d s 

Here, we assume that the top gate electric field is changed abruptly at the junction, such 
as E t = E tn for x < and E t = E tp for x > (abrupt junction). The & is constant 
in homogeneous regions which are sufficiently distant from x = 0. In these regions, the 
electric field in the semiconductor is zero (E- Y = —dfy/dx = 0), and the charge density 
is homogeneous (q(x) = e t E t — S\,E h ). A p-n junction is formed at x = for eN p = 
e t E tp — £t>-^b > and —eN n = e t E tp — e^E h < 0. There is a transient region form p to 
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n around x = 0. If the 2D conductor has the band gap, charge carriers are depleted from 
the transient region. Assuming that \q{x)\ is much smaller than eN p and eN n , the Poisson 
equation in the depletion region becomes, 



d 2 V eN p 
d 2 $ eN n 



dx 2 £ s gL 



(for the p side), 
(for the n side). 



These equations have same form as in a 3D p-n junctional], in which the acceptor density 
(Na) and the the donor density (N-d) are replaced by N p /d s and N n /d s , respectively. Use 
the same solution method for the 3D case, we obtain a depletion region width for the abrupt 
junction, 



abrupt 



2e s d s N p + N„ A g 



\ e N p N n 

where A $ is the electrical potential difference between the p side and the n side. Using 
typical values in our experiment (N p ~ N n ~ 2 x 10 12 cm -2 , e s = 3.0e for graphite[33j, 
d s = 0.67 nm for bilayer graphene, and ~ E g = 0.2 eV), the transient region width is 

^abrupt _ 2 nm _ 

The abrupt junction model is applicable only when the transient region, in which the 
top gate electric field changes from E tn to E tp , is narrower than the W^^. If not, we 
have to take the transient region width as the depletion region width H^dIS]. In the stepwise 
top gate, the top-gate dielectric thickness is changed from d% to d 2 {d\ < d 2 ) at position 
x = (thick solid line in Fig. A-2(a)). An electric potential, (f)(x,z), between graphene 
and the top gate is obtained by solving a Poisson equation numerically with a boundary 
condition, = at the graphene (z = 0) and = 1 at the top gate (z = d\ for x < 0, and 
z = g?2 for x > 0) (Fig. A-2(a)). In the calculation, SiCVequivalent thickness of the dielectric 
layer (d\ = 3.7 nm and d 2 = 6.3 nm) is used. The electric field component normal to the 
graphene, E± oc — d<ft/dz\ z= o, changes with a 5-nm-wide transient region around x = (Fig. 
A-2(b)). In the transient region, the electric field has a tangent component to the graphene, 
En oc —d(p/dx, at the vicinity of the graphene surface (Fig. A-2(c)). Since the transient 
region is wider than the PVp brupt (~ 2 nm), the transient region width gives the depletion 
region width, i.e., Wb ~ 5 nm. 
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Appendix B: Estimation of Band Diagram and Junction Resistance 

Equation (A-l) gives the charge density q(x) in graphene sandwiched by top and bottom 
gates, q(x) = e t E t — e^E h + e s d s dE-J dx. The gate electric fields are proportional the gate 
voltages, giving q(x) = —Ct g (x)Vt g —Cb g Vb g +e s d s dEi/dx+qo with C tg (bg) as local capacitance 
for the top (bottom) gate and go as the charge density for Vb g = V tg = 0. Three typical band 
diagrams for graphene with spatially modulated charge density are mentioned in Fig. B-l(a)- 
1(c), corresponding to p-p, n-p, and n-n combinations. In the resistance measurement, these 
combinations are separated by the two ridges of CNP (Fig. B-l(d)). On the CNPs, the 
Fermi level is at the mid-gap of the right or left side in the stepwise modulation. 

For the Vt g fixed in between the two CNPs, a spatial modulation forming a junction with 
p- and n-region is generated. When the Vbg is applied to the p-n junction, the p- and the 
n-region are separated by the tunneling barrier caused by the band gap (Fig. B-l(b)). The 
charge density in each p- (or n-) region is given by q p ( OTn ) = — C tg (V tg — VQ^ rn ^), where 
VcN° r is the top gate voltage for the CNP in the p(or n)-region. When the Vb g = —40 V 
and Vtg = 1-0 V (corresponding to "(b)" in Fig. B-l(d)) are applied, q p /e = 1.7 x 10 12 cm -2 
(or q n /e = 2.5 x 10 12 cm -2 ) is extracted from V^n = 1.5 V (or Vq N = 0.57 V). The carrier 
density q p (o rn )/e equals to integration of the density of states (DOS) from the mid-gap to 
the Fermi level. The typical charge density q p /e ~ q n /e ~2x 10 12 cm -2 in our experiment 
is similar to that in reported experiments on graphene p-n junction [3, 4J. 

An ideal DOS (D) without band tail is given as a function of the energy (e) measured from 
the mid-gap: D(e) ~ {t ± /c 2 )^A/\e - E g /2\ for |e| > E g /2 and D(e) = for |e| < E g /2, 
with A as the potential difference between the two graphene layers and c(~ 1 x 10 6 m/s) is 
the Dirac velocity in graphene [29]. The A is related to the potential difference between the 
two graphene layers by E g = 2(A — A 3 /tj_ 2 ), where tj_(~ 0.3 eV) is the interlayer hopping 
integraljSl], giving A ~ 0.12 eV for E g ~ 0.2 eV. Because the DOS of the graphene has a 
singularity at the band edges (e = ±E g /2) [2H1 12H], the Fermi level E p (or E n ) measured 
from the band edge of the p- (or n-) region is related to the charge density q p (or q n ) as 
qi/e = — jai(E s /2+Ei) jj^ds = — ai(2tj_/c 2 )y/A ■ Et (i = p or n), where a p = —1 and a n = 1. 
Using typical values of the carrier density of q p /e = q n /e ~ 2 x 10 12 cm -2 with potential 
difference between the two layers A ~ 0.12 eV and E g = 0.2 eV in our experiment, the 
Fermi level of the E p and E n in the each side of p-n junction are a few meV. These values 
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are two orders smaller than the energy gap E g , indicating that the E p and the E n are very 
close to the singularity peak in the each regions. In a realistic semiconductor graphene with 
band tails caused by localized states, the singularity peak becomes broader with increasing 
the density of localized states [28, 29J. Except for the graphene extremely disordered, the 
E p and E n are still much smaller than the E g , and are close to the band edges. Thus, the 
band diagram measured in the p-n junction can be described as Fig. B-l(b). 

The magnitude of the junction resistance can be roughly estimated from the band di- 
agram. The junction resistance consists of a parallel combination of a tunnel component 
and a leakage component. The tunnel component is estimated by the tunneling probabil- 
ity of an electron at the junction and the number of states involved in the tunneling. At 
the junction, the potential gradient is F ~ E g /Wn ~ 0.04 eV/nm, which leads to a tun- 
neling probability across the junction of T t = exp \-4y/2^E g 3/2 /(3hF)] ~ 0.07, where 
m* = t± 2 /4Ac 2 ~ 0.03m e for A = 0.12 eV, m e is the bare electron mass, and c ~ 1 x 10 6 
m/s is the Dirac velocity in graphene, for low energy states in a "Mexican hat" band 
structure [HI [331 EHl EH]- Landauers formula gives the tunnel conductance 1/Rt ~ N (e 2 / 7th)T t 
, where N is the number of quantized wavenumber states in the channel-width direction and 
h is Planck's constant. Even at the band edge, the Mexican hat band structure has a non- 
zero wavenumber range from to fco, where k$ = A/hc (~ 0.2 nm" 1 for A = 0.12 e V) [13| [54] . 
The range of the wavenumber in the p-type and n-type region is approximately to k , be- 
cause the Fermi level is close to the band edge. N is roughly k W/n ~ 30 (W ~ 0.4 /zm 
is the channel width) because the wavenumber is quantized [51)] by tt/W. Thus, the tunnel 
resistance Rt is estimated to be several kfl. The leakage conduction via localized states 
coexists with the tunnel conduction. The leakage resistance is expressed as i?L ~ Po^d, 
where po is the residual resistance at the mid-gap states; that is, the resistance per unit 
length at the CNP. Using p ~ 0.2 kfl/nm for V^ g = —40 V, the leakage resistance is on 
the order of ~ 1 kfl. Thus, the leakage current is comparable to the tunneling current. 
Finally, the junction resistance obtained, i?j = l/(l/i? t + 1/Rl), is on the order of 1 kfl. 
Using this value, it is estimated that the voltage drop at the junction Vj is a few percent of 
Vds, corresponding to a few mV, because the total resistance of the sample is about 40 kfl. 
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FIG. 1: (Color online) (a) Schematic diagram of p-n junction mechanism for a gate-controlled 2D 
semiconductor (left) and an impurity-doped 3D semiconductor (right), (b) Optical micrograph 
of a BLG p-n junction with a stepwise top gate. The solid lines show the shape of the graphene 
channel made by oxygen plasma etching. The dashed lines show the area in which the graphene 
is covered by a 5-nm-thick Si02 layer, (c) Schematic illustration of the cross-sectional view of the 
p-n junction. The locally inserted Si02 layer creates the stepwise structure in the top gate. The 
instrumental configuration for a four-terminal measurement is also shown. All transport properties 
shown in this paper were acquired in the same configuration, (d) Color map of the linear resistance, 
i?o, as a function of Vb g and V tg at T = 80 K, determined at ~ 1 mV. The dashed lines show 
the ridges of the CNPs separating p and n regions, (e) Schematic view of the geometrical structure 
of the stepwise top gate (top), the profile of the top gate electric field (middle), and the charge 
distribution (bottom). The top gate electric field is estimated numerically (Appendix A). 
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FIG. 2: (Color online) (a) Differential resistance, dV^s/dl , as a function of the source-drain 
voltage Vds for BLG under a relatively large electric field (Vbg = —40 V). The arrow indicates the 
dVds/dl peak. The top left inset shows current- volt age curve corresponding to the main panel. 
The bottom-right inset shows a dVd s /dI-Vd s curve for a relatively small electric field(Vbg = — 32 
V). (b) Distribution of the Vd s values for the dV^/dl peaks as a function of Vbg with various values 
of Vt g . The red closed symbols correspond to the data shown in (c). (c) dVd s /dI as a function of 
Vds and V"b g . The V g value is selected along the middle point of two ridges of CNPs in Fig. 1(c). 
All data were acquired at T = 80 K. 
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FIG. 3: (a) Schematic diagram of the DOS for the p-type and n-type sides: i) for the unbiased 
condition (Vj = 0) and ii)— iii) for the forward-biased (Vj > 0) conditions. The gray shaded areas 
show the density of occupied states with a Fermi distribution. Electron tunneling takes place from 
occupied states on one side to unoccupied states on the other side, (b) Schematic I-Vd s curve for 
the tunneling current (dashed curve) , the current via localized states (dotted curve) , and the total 
current (solid curve). The three conditions in (a) ( i), ii), and iii) ) are identified on the tunneling 
current curve and the total current curve. 
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FIG. 4: (Color online) (a) Temperature dependence of the dV& s /dI curve for Vb g = — 40 V and 
Vt g = 1.0 V. (b) Schematic diagram of A(E) = D n (E)D p (E + eVj) under a forward voltage 
drop of Vj corresponding to the tunneling current peak (Fig. 3(a)ii)). (c) Schematic current- 
voltage characteristics for the tunneling current at low temperature (thick dashed curve) and high 
temperature (thin dashed curve) and the total current with an additional leakage current (Fig. 3(b)) 
at low temperature (thick solid curve) and high temperature (thin solid curve), (d) Arrhenius plot 
of the dVds/dl values extracted from a at fixed values of V"d s for a forward bias (Vd s > 0). The 
dVd s /dI values are normalized at T = 80 K. The solid lines represent fittings within second order 
of the temperature. The dV^ s /dI curve is shown for a forward bias at T = 80 K in the inset of (b), 
in which the markers indicate the corresponding Vds values in the main panel. 
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FIG. A.l: Schematic illustrations of (a) 2D semiconductor sandwiched between a top gate and a 
bottom gate and (b) electric field in the semiconductor plane (E[) and in gate insulators (Et and 
E h ) 
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FIG. A. 2: (a) Contour plot of electric potential <fi with a boundary condition that <j> = at the 
graphene and = 1 at the stepwise top gate. The contour interval is 0.1. (b) Normal component 
of the electric field at the graphene surface, (c) Tangent component of the electric field at the 
vicinity of the graphene surface (z = 0.1 nm). 
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FIG. B.3: (Color online) Schematic band diagrams for combinations of (a) p-p, (b) n-p, and (c) 
n-n in the semiconductor graphene with the stepwise top gate, (d) The left panel shows resistance 
change as a function of the top gate voltage at fixed Vbg > an d the right panel shows a 3D plot 
of the resistance as a function of the bottom gate and top gate voltage for the same data in the 
middle panel of Fig. 1(c). The left panel corresponds to the cross section in the right panel. The 
marked three points in the right panel correspond to the band diagrams (a)-(c). 



